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ABSTRACT 

Aims. To account for the observed differential metallicity distribution (DMD) of the Milky Way halo, a semi-analytical model 
is presented in the framework of the hierarchical merging paradigm for structure formation. 

Methods. It is assumed that the Milky Way halo is composed of a number of sub-haloes with properties either as observed in 
the dwarf satellite galaxies of the Local group (shape of metallicity distribution, effective yield) or derived from calculations of 
structure formation (sub-halo distribution function). 

Results. With reasonable assumptions for the parameters involved, we find that the overall shape and effective yield of the 
Galactic halo DMD can be reproduced in the framework of such a simple model. The low metallicity tail of the DMD presents a 
defRciency of stars with respect to the simple model predictions (akin to the G-dwarf problem in the solar neighborhood); it is 
suggested that an early infall phase can account for that problem, as well as for the observed DMDs of dwarf satellite galaxies. 
Conclusions. Accretion of galaxies similar (but not identical) to the progenitors of present day dwarf satellites of the Milky Way 
may well have formed the Galactic halo. 

Key words. Galaxy: abundances, halo, formation; Galaxies: general, formation. Local Group 



1. Introduction 

' The metallicity distribution of long-lived stars is one of 
] the most powerful probes of galactic chemical evolution. 
' Indeed, the shape of the metallicity distribution is inde- 
' pendent of the rate at which stars are formed. It depends 
essentially on the stellar initial mass function (IMF), 
which fixes the "intrinsic" (true) yield of the stellar pop- 
ulation and the loss or gain of (gaseous) mass from the 
galaxian system (outflow or inflow, see next section and 
Prantzos 2007a). In the case of the solar neighborhood, 
for instance, the shape of the metallicity distribution sug- 
gests that the local disk was constructed by infall on long 
timcscales (e.g. Boissier and Prantzos 1999). 

The regular shape of the metallicity distribution of 
the Milky Way (MW) halo can readily be explained by 
the simple model of galactic chemical evolution with out- 
flow, as suggested by Hartwick (1976, see also Sec. 2.1). 
However, that explanation lies within the framework of 
the monolithic collapse scenario for the formation of the 
MW (Eggen, Lynden-Bell and Sandage, 1962). It seems 
dificult to keep that interpretation in the currently domi- 
nant paradigm of hierarchical galaxy formation. A coher- 
ent interpretation should be sought in the framework of 
the new paradigm, according to which the MW halo was 



formed from the hierarchical merging of smaller compo- 
nents (sub- haloes). 

Several attempts to account for the metallicity distri- 
bution of the MW halo in the modern framework were un- 
dertaken in recent years, starting with the work of Bekki 
and Chiba (2001). They model a galaxy closely ressem- 
bling the MW and despite the low spatial resolution of 
their simulations (2 kpc), they are able to provide a halo 
metallicity distribution with a peak at the observed metal- 
licity ([Fe/II]=-1.6), but with a shape which is much nar- 
rower than the observed one ; theyobtain a double peaked 
metallicity distribution (which is not observed), proba- 
bly reflecting late major mergers due to the small mass 
resolution. In their five MW-type simulations (in a cos- 
mological framework and using various prescriptions for 
onflow/feedback) Scanapieco and Broadhurst (2001) find 
no ressemblance of the resulting halo metallicity distribu- 
tions with the observed one. This is also the case in Font et 
al. (2006), who find systematically too metal rich haloes 
with respect to the MW halo (see their Fig. 6). On the 
contrary, using semi-analytical models (starting with the 
extended Press-Schecter formalism and using various pre- 
scriptions for star formation, feedback and mixing of the 
stellar ejecta) Tumlinson (2006) and Salvadori, Schneider 



2 



N. Prantzos: HMD 



and Ferrara (2007) find rather good agreement between 
the model and observed halo metallicity distribution. 

Independently of their success or failure in reproduc- 
ing the observations, those recent models provide little or 
no physical insight into the physical processes that shaped 
the metallicity distribution of the MW halo. Indeed, if the 
halo was built from a large number of sucessive mergers 
of sub-haloes, why is its metallicity distribution so well 
described by the simple model with outflow (which refers 
to a single system)? And what determines the peak of the 
metallicity distribution at [Fe/H]='^-1.6; which is (suc- 
cessfully) interpreted in the simple model by a single pa- 
rameter (the outflow rate) ? Even in the succesfuU models 
cited in the previous paragraph, no clear answer is pro- 
vided to these questions. 

In this work we present an attempt to built the halo 
metallicity distribution analytically, in the framework of 
the hierarchical merging paradigm. Preliminary results 
have already been presented in Prantzos (2007b). It is as- 
sumed that the building blocks were galaxies with proper- 
ties similar to (but not identical with) those of the Local 
group dwarf galaxies that we observe today. In that way, 
the physics of the whole process become (hopefully) more 
clear. The plan of the paper is as follows: In Sec. 2 we 
present the simple model of the outflow and its overall 
successfull description of observations; we discuss, in par- 
ticular, some issues related to its low- and high- metallicity 
tails. In Sec. 3 we present the properties of the sub-haloes 
that formed the MW halo. These are derived either from 
observations of the Local froup dwarf galaxies (the shape 
of the metallicity distribution. Sec. 3.1 and the effective 
yield. Sec. 3.2) or from the theory of large scale struc- 
ture formation (the distribution function of baryonic sub- 
haloes. Sec. 3.3). The results of the model are presented 
in Sec. 3.4 and summarized in Sec. 4. 

2. The halo metallicity distribution and the simple 
model 

The halo metallicity distribution for field stars was well 
established in the metallicity range -2.2 < [Fe/H] < - 
0.8 by two major surveys in the 90ies (Ryan and Norris 
1991, Carney et al. 1996). An ongoing survey, known as 
the Hamburg/ESO survey, seeks to establish the form of 
the metallicity distribution at even lower metallicities. Its 
results should be considered as representative down to 
[Fe/H] '^-3, and thus they superseed the survey of Ryan 
and Norris (1991) below [Fe/H]=-2. For [Fe/H]>-2.2, var- 
ious biases affect the results of the Hamburg/ESO survey. 
Fig.l displays a composite of the Ryan and Norris (1991) 
and the Hamburg/ESO surveys, adjusted to fit each other 
at [Fe/H] =-2. 2. In this work we shall consider the main 
body of that distribution, and we shall only bricfiy discuss 
some issues related to its low- and high-metallicity tails. 

As already noticed a long time ago (Hartwick 1976) 
the halo metallicity distribution is nicely described by the 
simple model of galactic chemical evolution (GCE). In the 



framework of that model, the metallicity Z is given as a 
function of the gas fraction [i as 

Z ^ Vln(^ + ^0 (1) 

where Zq is the initial metallicity of the system (zero for all 
isotopes, except those created in the Big Bang) and p is the 
yield of a stellar generation, i.e. the newly created mass of 
a metal returned to the ISM per unit mass blocked in stars 
(both long lived stars and stellar remnants). Metallicities 
and yield are expressed in units of the solar metallicity Zq. 
The yield depends on the stellar nucleosynthesis models 
and on the adopted stellar initial mass function (IMF). If 
the system evolves at a constant mass (closed box), the 
yield is called the true yield, otherwise (i.e. in case of mass 
loss or gain) it is called the effective yield] its value is then 
always lower than the one of the true yield and it can 
be analytically evaluated in some specific cases (see e.g. 
Edmunds 1990). 

Taking into account that the number of stars is pro- 
portional to the stellar mass (for a constant IMF) and 
that the star mass fraction is given by — 1. — /x, 
one may use Eq. (1) to derive the differential metallicity 
distribution (DMD) of the system as: 



d{n/ni) 
dilogZ) 




where Zi is the final metallicity of the system and rii the 
total number of stars (having metallicities < Zf). This 
function has a maximum for Z — = p, allowing one to 
evaluate easily the effective yield p if the DMD is deter- 
mined observationally. 

The DMD of field halo stars peaks at [Fe/H]^ — 1.6, sug- 
gesting an effective halo yield p^aio ~l/40 Zq for Fe. This 
has to be compared to the true yield obtained in the so- 
lar neighborhood. The peak of the local DMD is located 
in the range -0.2<[Fe/H]<0., depending on the adopted 
metallicity calibration (see e.g. Haywood 2006, Holmberg 
et al. 2007). Assuming that the peak is at [Fe/H] =-0.1, it 
points to a local yield of p£)isfc=0.8 Zq. This value is not 
very different from the yield derived for the bulge of the 
MW, which evolved more or less as a closed box (as sug- 
gested from its metallicity disribution, see e.g. Fulbright 
et al. 2006); this is considered to be the true yield (as- 
suming the same IMF between local disk and the bulge). 
Taking into account that ~60-65% of solar Fc are pro- 
duced by SNIa and only ~35-40% by SNII (e.g. Goswami 
and Prantzos 2000) it turns out that the true Fe yield 
of SNII during the local disk evolution was poisk ~0.32 
Zq and, consequently, the effective Fe yield of SNII during 
the halo phase (where they dominated Fe production) was 

PHalo 

'^O.OS PDisk- 

The simplest interpretation of such a difference be- 
tween the effective yields in the halo and the local disk 
remains still the one of Hartwick (1976), who suggested 
that outflow at a rate F ^ k (where ^ is the Star 
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Formation Rate or SFR) occurcd during the halo forma- 
tion. In the framework of the Simple model of GCE such 
outflow reduces the true yield pxrue to its effective value 
PHaio = PTrue{l-R)/{l + k-R) (e.g. Pagcl 1997) where R 
is the Returned Mass Fraction [R ~0. 32-0. 36 for most of 
the modern IMFs, like e.g. the one of Kroupa et al. (1993) 
or of Kroupa 2002). The halo data suggest then that the 
ratio of the outflow to the star formation rate was 



k = {l-R) 



PTr 



- 1 



PHalo 



(3) 



and, by replacing pTrue with poisk one finds that k ~7-8 
times the star formation rate. This high outflow rate could 
be interpreted cither as gas expulsion after heating by su- 
pernova explosions, or as simple flowing of gas through the 
system, e.g. towards the Galaxy's bulge. 

Although the siple model of GCE with outflow fits ex- 
tremely well the bulk of the halo DMD, the situation is less 
satisfactory for the low metallicity tail of that DMD. The 
problem (a lack of low metallicity stars w.r.t. the predic- 
tions of the simple model) was noticed sometime ago and 
prompted some interesting scenarii concerning the physics 
of the early galaxy, like e.g. pre-enrichment by a puta- 
tive population of Pop. Ill stars (Norris 1999) or various 
possibilities of inhomogeneous early evolution (Tsujimoto 
et al. 1999, Oey 2003). Prantzos (2003) showed that the 
situation is in fact even worse if the instantaneous recy- 
cling approximation (adopted in all previous works on the 
problem) is dropped and the finite lifetimes of stars are 
taken into account. Noting that the situation is reminis- 
cent of the G-dwarf problem in the solar neighborhood, he 
further suggested that a similar solution should naturally 
apply, namely an early phase of rapid infall (in a time 
scale of less than 0.1 Gyr) forming the Milky Way's halo. 
Such a scenario is illustrated in Fig. 1. It may account for 
the "G-dwarf problem of the halo, at the price of intro- 
ducing more degrees of freedom than the simple model. 
At present, the situation is not quite clear regarding the 
low metallicity tail of the halo DMD. We shall present 
some tentative conclusions, based on preliminary data, in 
Sec. 3.5; however, more sound conclusions should be drawn 
only after the definitive results of the full Hamburg/ESo 
survey (Schork et al. in preparation). 

The high metallicity tail of the halo DMD may also 
hold some important clues as to the halo evolution. Its 
shape is surprisingly well fitted by the simple model, up 
to [Fe/H]~-1. In the framework of the simple model, this 
implies that the halo evolved down to a very low gas 
fraction /^(< Zi) ~0.05. Indeed, since log(Zi)~-l and 
log(pi/a;o)~"l-6, the main result of the simple model (Eq. 
1) leads to ji ~0.05. However, a different interpretation 
was suggested by Oey (2003). In the framework of her 
Simple Inhomogeneous Model, the halo is a relatively un- 
evolved system, composed of a number of regions enriched 
to different levels (the larger ones having lower mctallici- 
ties). In that model, the final metallicity distribution re- 
flects essentially the size distribution of those regions, and 
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Fig. 1. Differential metallicity distribution of field halo 
stars in linear {top) and logarithmic {bottom) scales. Data 
are from Ryan and Norris (1991, open symbols)^ and the 
ongoing Hamburg/ESO project ( filled symbols). The two 
data sets are normalised at [Fe/H]-2.2; above that value, 
the Hamburg/ESO data are incomplete, while below that 
metallicity the Ryan and Norris (1991) data set is incom- 
plete. The dotted curve is a simple model with instanta- 
neous recycling (IRA) and outflow rate equal to 8 times 
the star formation rate. The solid curve is obtained as 
in Prantzos (2003), from a model without IRA, an early 
phase of rapid infall and a constant outflow rate equal to 
7 times the SFR. All curves and data are normalised to 
max= 1 . 



at high metallicitics it scales as Z ^, in satisfactory agree- 
ment with observations. 

At present, it is hard to distinguish between the two 
alternatives, i.e. to decide whether the strong decrease of 
the observed star number in the high metallicity tail of the 
halo is due to the fact that the system is highly evolved 
or to its inhomogeneous nature. Obviously, it cannot be 
due to observational errors alone, since the width of the 
DMD (~1 dex) is much larger than the one justified by 
such errors(0. 1-0.2 dex). We simply note that in the mod- 
ern framework of galaxy formation through hierarchical 
merging, smaller units have smaller metallicitics, as dis- 
cussed in the next section. 
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3. The halo DMD and hierarchical merging 

In this work we assume that the MW halo was formed 
by the merging of smaller units ("sub-haloes"), as implied 
by the hierarchical merging scenario for galaxy formation. 
It should be noticed at this point that, despite its over- 
all sucess with large scale structure in the Universe, this 
scenario still faces important problems when it comes to 
predictions concerning the behaviour of baryons at galac- 
tic scales. Among those problems, the one concerning the 
formation chronology of baryonic vs dark matter struc- 
tures (namely the fact that large galaxies are observed to 
be older than small ones, while large dark matter haloes 
assemble later than low mass ones) is known as down- 
sizing. The word downsizing was first used by Cowie et 
al. (1996) on observational grounds, while the problem 
with spirals was noticed in Boissier and Prantzos (2000) 
and for ellipticals in Thomas et al. (2002). van den Bosch 
(2002) showed explicitly how semi-analytical models in 
the framework of hierarchical merging scenarios fail to re- 
produce the colour-magnitude relation of galaxies, even 
if some form of feed-back is included. Despite several at- 
tempts, no universally acepted model meeting the obser- 
vational requirements has been proposed up to now (but 
see Cattaneo et al. 2008, for a possible solution in the case 
of elliptical galaxies). 

Despite this important shortcoming, several models 
have been built in the framework of the hierarchical merg- 
ing scenario, exploring its consequences on the formation 
of the MW halo and some of them address specifically the 
issue of the halo metallicity distribution (e.g. Tumlinson 
2006, Salvadori et al. 2007). In these models star formation 
takes place early on in the different baryonic sub- haloes, 
which subsequently are tidally disrupted to form a halo 
component with an old stellar population. Such an early 
star formation is required to account for the observed age 
of halo stars (>11 Gyr, e.g. Hill et al. 2002) and their large 
[Q;/Fe]~0.5 ratio (Cayrel et al. 2004) , suggesting that those 
stars had not been enriched by Fe from SNIa. 

We shall assume then that the MW sub-haloes formed 
their stars early on, plausibly because at those early times 
they evolved already in a dense environment (i.e. the deep 
potential well of the MW halo); similar size sub-haloes 
in less dense environments should form their stars at a 
slower pace, on average, and develop thus young stellar 
populations today. During the short period of star foma- 
tion in each of the sub-haloes, their shallow individual 
potential wells could hardly keep their gas (affected by 
ram-pressure stripping and heated by supernovae) which 
would then leave the sub-haloes (probably to end in the 
bulge, as suggested long ago on the basis of angular mo- 
mentum conservation arguments, e.g. Gilmore et al. 1989). 
Obviously, the amount of outflowing gas (and of the met- 
als kept inside the system, i.e. the effective yield) should 
depend on the depth of the potential well of each sub- 

^ In fact, over the years, the word downsizing was used to 
denote several diflFerent, albeit not quite unrelated, things (see 
Conroy and Weschler 2008). 



halo, i.e. on the mass of the corresponding dark matter 
sub-halo. 

In order to calculate semi-analytically the resulting 
DMD in the framework of the aforcmcntionncd scenario, it 
is further assumed that each of the sub-halocs had a DMD 
described well by the simple model, i.e. by Eq. (2). It re- 
mains then to evaluate the corresponding effective yield 
p{M) of each sub-halo as a function of the mass of the 
sub-halo, as well as the mass function of the sub-haloes 
dN/dM, in order to obtain the required result. Each one 
of those ingredients is discussed in the following sections. 

3.1. The shape of sub-halo DMD 

It is assumed here that each one of the merging sub-haloes 
has a DMD described by the simple model with an appro- 
priate effective yield. This assumption is based on recent 
observations of the dwarf spheroidal (dSph) satellites of 
the Milky Way, as will be discussed below. It is true that 
the dSphs that we see today cannot be the components of 
the MW halo, because of their observed abundance pat- 
terns: indeed, as discussed in several places (e.g. Shetrone, 
Cote and Sargent 2001; Venn et al. 2004) their a/Fe ra- 
tios are typically smaller than the [a/Fe]^0.4^const. ra- 
tio of halo stars. This implies that they evolved on longer 
timescalcs than the Galactic halo, allowing SNIa to enrich 
their ISM with Fc-pcak nuclei and thus to lower the a/Fc 
ratio by a factor of ~2 (as evidenced from the [O /Fc]~0 ra- 
tio in their highest metallicity stars, about half of the MW 
halo value). However, the shape of the DMD of the simple 
model does not depend on the star formation history or 
the evolutionary timescale. This is precisely what makes 
the DMD such a powerful diagnostic tool: it provides ac- 
curate and indepedent information on gas flows into and 
out of the system, as well as on the initial metallicity. In 
other terms, the question is whether the components of 
the MW halo had a DMD dccribcd by the simple model 
or by some different form. To answer that question, one 
may turn cither to simulations (which arc difficult to con- 
strain) or to observations of galaxian systems smaller than 
the MW halo. This is the case of the Local group dSphs. 

Helmi et al. (2006) report observations of several hun- 
dred red giant stars in four nearby dSphs (Sculptor, 
Fornax, Sextans and Carina) , obtained with the FLAMES 
spectrograph in the VLT. The surface [Fe/H] ratio of those 
stars should only marginally be affected by the 1st dredge- 
up, and thus they may be used as tracers of the chemi- 
cal evolution of the corresponding galaxies. The resulting 
DMDs are displayed as histograms in Fig. 2, where they 
are compared to the simple model with appropriate effec- 
tive yields {solid curves). The effective yield in each case 
was simply assumed to equal the peak metallicity (Eq. 2). 
It can be seen that the overall shape of the DMDs is quite 
well fitted by the simple models. This is important, since 
i) it strongly suggests that all DMDs of small galaxian 
systems can be described by the simple model (at least, 
there are no counter examples to that hypothesis) and ii) 
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Fig. 2. Metallicity distributions of dwarf satellites of the 
Milky Way. Data are in histograms (from Helmi et al. 
2006). Solid curves indicate the results of simple GCE 
models with outflow proportional to the star formation 
rate; the corresponding effective yields (in Zq) appear on 
top right of each panel. Dashed curves are fits obtained 
with an early infall phase, while dotted curves are models 
with an initial metallicity log(2'o) ^-3; both modifications 
to the simple model (i.e. infall and initial metallicity) im- 
prove the fits to the data. 

it allows to determine effective yields by simply taking the 
peak metallicity of each DMD ; this will be used in Sec. 
3.2 below. 

Before proceeding to the determination of effective 
yields, we note that the fit of the simple model to the 
data of dSphs fails in the low metallicity tails: as already 
noticed in Helmi et al. (2006), there is a lack of low metal- 
licity stars w.r.t. the simple model predictions. Helmi et al. 
(2006) attribute this to a pre-enrichment of the gas out of 
which the dSphs were formed, i.e. they assume that Zq ^0 
in Eq. (1) and (2). By fitting their data that way, they de- 
termine then initial metallicities in the range [Fe/H]o ~- 
2.9 - -2.7 for the four galaxies. Furthermore, they argue 
that this result has a strong implication ; it would im- 
ply that "... the progenitors of the Milky way (halo) and 
of dSphs must have been different" . However, this is not 
necessarily true. 

Indeed, the idea of a finite initial metallicity (pre- 
enrichment) was for a long time considered as the "sim- 
plest" solution to the G-dwarf problem in the solar neigh- 
borhood (see Pagel 1997). But, although mathematically 
the simplest, it was not the most "natural", and infall has 
now superseed it. In the case of the MW halo, the situation 
is less clear, but early infall appears again as an interest- 
ing scenario, as argued in Prantzos (2003) and Sec. 3.1. In 



the light of these arguments, one should not exclude that 
early infall is also the solution to the " low metallicity tail" 
problem of the DMDs of the local dSphs. If an early infall 
phase turns out indeed to be a universal feature of galax- 
ian systems, then the MW halo may well be composed of 
galaxies similar to the progenitors of the dSphs. We shall 
return to the low metallicity tail of the halo in Sec. 3.5, 
after considering first the bulk of the halo DMD. 

3.2. The effective yield of the sub-haloes 

If the DMDs of all the components of the Galactic halo 
are described by the simple model, then their shape is es- 
sentially described by the corresponding effective yield p 
(and, to a lesser degree, by the corresponding initial metal- 
licity Zq). Observations suggest that the effective yield is 
a monotonically increasing function of the galaxy's stellar 
mass M*. Data for about 40 galaxies of the Local group, 
compiled by Dekel and Woo (2003) are presented in Fig. 
3. It can be seen that the (median) stellar metallicity of 
those galaxies varies as 0.01 (M^/IO^Mq)"-''. Dekel and 
Woo (2003) offer a plausible theoretical framework for that 
corelation, based on supernova feedback and larger mass 
loss from smaller galaxies; their model explains satisfacto- 
rily other observed properties, as e.g. velocity dispersion 
and central surface brightness vs. stellar mass. This is also 
the relationship between the effective yield (equal to the 
median stellar metallicity from Eq. 2) and the mass of 
those galaxies. In the case of the progenitor systems of 
the MW halo, however, the effective yield must have been 
lower, since SNIa had not time to contribute (as evidenced 
by the high a/Fe~0.4 ratios of halo stars), by a factor of 
about two (see discussion in previous section) . We assume 
then that the effective yield of the MW halo components 
(accreted satellites or sub-haloes) is given (in solar units 
Zq) by 



0.005 



0.4 



(4) 



i.e. the thick dotted curve in Fig. 3. 

Obviously, the stellar mass M^, of each of the sub- 
haloes should be < Mh where Mh is the stellar mass 
of the MW halo. In a recent paper. Bell et al. (2007), 
analysing ~4 10^ stars of the MW halo find that the stel- 
lar mass in the halo region at distances 1 to 40 kpc from 
the Galactic center is Mh= 4±0.8 10^ Mq and we adopt 
that value here. The conclusions of our work depend only 
slightly on the precise value of Mh. 

An inspection of Fig. 3 shows that the effective yield of 
the MW halo is substantially lower (factor '^4) than the 
one of dwarf galaxies with stellar masses in the range of 
Mh- Part of this difference (factor ~2) is due to the fact 
that no SNIa contributed to the effective Fe yield of the 
MW halo. The remaining difference can be qualitatively 
understood if it is assumed that the halo was made by 
the assembly of smaller mass satellites: the resulting total 
mass {Mh) is larger than the mass of the largest sub- 
halo, but the resulting mass-weighted effective yield {ph) 
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Fig. 3. Stellar metallicity vs stellar mass for nearby galax- 
ies; data and model {upper curves) are from Dekel and 
Woo (2003), with dl standing for dwarf irregulars and 
dE for dwarf ellipticals. The thick dotted line represents 
the effective yield of the sub-haloes that formed the MW 
halo according to this work (i.e. with no contribution from 
SNIa, see Sec. 3.2). The MW halo, with average metallicity 
[Fe/H]=--1.6 and estimated mass ~4xl0* Mq falls below 
both curves. 



is smaller than the effective yield of the largest sub-halo. 
On the contrary, if the MW halo were formed as a single 
system, it is hard to understand why its effective yield is 
so low w.r.t. systems of similar mass. 

3.3. The stellar mass function of the sub-haloes 

Contrary to the cases of the effective yield and of the shape 
of the DMD, the mass function of the sub- haloes which 
presumably formed the MW halo can not be observed to- 
day (since the observed number of local dwarf galaxies is 
insufficient for that purpose) but has to be derived from 
semi-analytical or numerical calculations of structure for- 
mation. 

Hierarchical galaxy formation scenarios predict the 
mass function of the dark matter sub-haloes which com- 
pose a dark matter halo at a given rcdshift. In the high 
resolution simulation Via Lactea, one of the most "real- 
istic" simulations of a MW-like galaxy performed so far, 
Diemand et al. (2007) find that the cumulative dark sub- 
halo mass function at redshift z=0 can well be approxi- 
mated by a power-law iV(> M£,)=0.0064 {Md/Mdh)~'^ 
(Mj) being the mass of the dark matter sub-haloes and 
MjjH the total mass of the dark halo); this approxima- 
tion is valid over three orders of magnitude, i.e. from 3 
10^ to 3 10^ M0 inside the galactic virial radius. 



The corresponding differential mass function is obvi- 
ously 



dM^ - 



(5) 



Salvadori et al. (2007) find in their numerical simulations 
that the same result is valid at higher redshifts (even up to 
z=5, see their Fig. 3), at least for the low mass part of the 
dark halo mass function. In a recent work, Giokoli et al. 
(2007) find analytically (by applying the extended Press 
Schectcr formalism and using the progenitor (conditional) 
mass function) that the result of Eq. (5) is independant 
of the mass of the progenitor halo and extends down to 
sub-halo masses as small as 10~^ of that mass. 

However, in our case, we are interested in the mass 
function of the stellar sub-haloes, and not of the dark 
ones. The initial ratio of baryonic mass Mg to dark mat- 
ter mass is usually taken to be independent of the dark 
halo mass Mo and equal to the corresponding cosmic ra- 
tio f = VLb/^d ~0.17, where VLb '-0.04 and Vto ~0.23 
are the cosmic densities of baryons and dark matter, re- 
spectively (e.g. Fukugita and Peebles 2003). If that ra- 
tio remained constant throughout galaxy evolution, then 
Eq. (5) would obviously provide the required stellar mass 
function also (assuming that all, or most of, baryons form 
stars). However, several effects may subsequently alter 
that ratio, like e.g. ram pressure stripping or galactic 
winds from supernova feedback. Observations and theoret- 
ical arguments suggest that those effects become more im- 
portant in lower mass systems (e.g. Dekel and Woo 2003). 
It is difficult to evaluate the impact of those effects on 
the final baryonic fraction of a given galaxy, and to derive 
the corresponding mass function of baryonic (mostly stel- 
lar) sub-haloes just from first principles. We attempt here 
an analytical evaluation of those quantities, based on the 
various empirical relations established in Sec. 2 and 3. 

From Eq. (3), one can see that the ratio of the outflow 
to the star formation rate from a galaxy (or a sub-halo) 
of stellar mass Af* is 



(6) 



since pTrue >> p{M^) for those small galaxies/sub-haloes, 
which are dominated by outflows (as can be seen in Fig. 4, 
top panel, even for a galaxy of ~10* Mq, the outflow 
rate is 3 times the SFR). By virtue of Eq. (4), one has 
then fc(M*) oc M"" **, in other terms 



M, 



9EL oc M-°-4 



(7) 



where Mqut = / Fdt ™ / k'^dt is the mass lost from the 
system and Af* = / ^'di is the stellar mass of the system. 
On the other hand, the barionic mass Mb is given by 



MouT + M^ = Mb = / Md 
with / = CIb/^d- 



(8) 
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Fig. 4. Properties of the sub-haloes as a function of tlicir 
stellar mass, empirically derived as discussed in Sec. 3. 
From top to bottom: Outflow rate, in units of the cor- 
responding star formation rate; Effective yield, in solar 
units; Distribution function; Cumulative fraction of stel- 
lar mass contributed by the sub- haloes. The total mass of 
the MW halo is 4 10® Mq. 

For MouT » M*, Eqs (7) and (8) lead to M°-^ oc 
Md ■ Combining this with the distribution function of the 
dark matter sub-haloes (Eq. 5), one finally gets 



dN 



oc m: 



(9) 



i.e. the distribution function of the stellar sub-haloes is 
flatter than the distribution function of the dark mat- 
ter sub-haloes because of the mass loss, which becomes 
progressively stronger as we move to less massive sub- 
haloes. The normalisation of the stellar sub-halo distri- 
bution function is made through 



Ah 



Ml 



dN 
dAL 



H 



(10) 



where M//=4 10® M0 is the adopted stellar mass of the 
MW halo. 

The lower mass limit Mi is adopted here to be Mi =1-2 
10^ Mq, in agreement with the lower mass bound of ob- 
served dSphs in the Local group (see Fig. 3). Such galax- 
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Fig. 5. Top and middle panels: Differential metallicity dis- 
tribution (in lin and log scales, respectively) of the MW 
halo, assumed to be composed of a population of smaller 
units (sub-haloes). The individual DMDs of a few sub- 
haloes, from 10^ Mq to 4 lO"^ Mq, are indicated in the 
middle panel, as well as the sum over all haloes {solid 
upper curves in both panels, compared to observations). 
Dotted curves in top and middle panels indicate the re- 
sults of the simple model with outflow (same as in Fig. 1). 
Because of their large number, small sub-haloes with low 
effective yields contribute the largest fraction of the lowest 
metallicity stars, while large haloes contribute most of the 
high metallicity stars {bottom panel). 



ics have internal velocities V >10 km/sH Dekel and Woo 
(2003) argue that the gas in haloes with V <10 km/s can- 
not cool to form stars at any early epoch and that dwarf 
ellipticals form in haloes with 10< V <30 km/s. From 
their V vs M* diagrams (their Figs. 3 and 5) one sees 
that a velocity of 30 km/s corresponds to a stellar mass 
M» ^2 10® Mq, i.e. about half the adopted stellar mass of 
the MW halo. This is the highest conceivable mass for a 
component of a ~4 10® Mohalo and it is the upper mass 
limit M2 that we adopt here for Eq. (10). 



^ Dekel and Woo (2003) calculate internal velocity as V 
\/3ap, where ap is the projected central velocity dispersion. 
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The main properties of the sub-halo set constructed in 
this section appear in Fig. 4 as a function of the stellar 
sub-halo mass M^. 

3.4. The composite stellar halo DMD 

We now have all the ingredients required to calculate the 
DMD of the MW halo, assuming that 

1) It is composed of a number of sub- haloes with stellar 
masses ranging from A/i=2 10^ Mq to M2=2 10® Mq and 
distributed as in Eq. (9); and 

2) Each of the stellar sub-haloes has a DMD given by 
the simple model (Eq. 2) with a mass dependent effective 
yield given by Eq. (4). 

The resulting DMD is then obtained as a sum over all 
sub-haloes: 



d{n/ni) 
d{logZ) 



AU 



Ml 



d[n(M,)/ni(M,)] dN 



dilogZ) 



dAL 



M, dM, (11) 



The result appears in Fig. 5 (top panel in linear and middle 
panel in logarithmic scales, respectively). It can be seen 
that it fits the observed DMDs at least as well as the 
simple model a la Hartwick. In fact, at low metallicities, 
the behaviour of the simple model (with outflow rate = 
8 times the SFR) is similar to the one of the composite 
model. The reason is that far from the peak (given by the 
effective yield) the DMDs of all individual haloes have the 
same slope (see middle panel of Fig. 5) and so does their 
sum. In the high metallicity region, the composite model 
does even a little better than the simple model. This region 
is dominated by the DMDs of the most massive haloes, 
which peak at [Fe/Hj^-l, whereas at low metallicities the 
small haloes dominate the composite DMD (see bottom 
panel in Fig. 5). 

In summary, under the assumptions made in this sec- 
tion, the bulk of the DMD of the MW halo results nat- 
urally as the sum of the DMDs of the component sub- 
haloes. It should be noted that all the ingredients of the 
analytical model are taken from observations of local satel- 
lite galaxies, except for the adopted mass function of the 
sub-haloes (which results from analytical theory of struc- 
ture formation plus a small modification to account for the 
role of outflows). Obviously, by assuming different values 
for the slope of the dark matter haloes in Eq. (5) and for 
the mass limits Mi and A/2 in Eqs. (10) and (11), one may 
modify the peak of the resulting composite DMD, thus al- 
lowing for differences between the halo DMDs of different 
galaxies. 

3.5. The low metallicity tail of the DMD 

As already discussed in Sec. 2, current surveys of the low 
metallicity tail of the halo DMD flnd a lack of stars w.r.t. 
the predictions of the simple model (with outflow) be- 
low [Fe/H]'^-3. A similar trend is found in satellite galax- 
ies of the Milky way by Helmi et al. (2006), albeit at 
slightly higher metallicities ([Fe/H]~-2.8). In the case of 



the MW halo, Prantzos (2003) argued that an early in- 
fall phase may account for the discrepancy (see Fig. 1). 
Obviously, the same solution may also apply to the case of 
the MW satellites, i.e. the conclusion of Helmi et al. (2006) 
that those galaxies were formed from pre-enriched gas (of 
[Fe/II]~-2.8) does not necessarily hold. Consequently, the 
main point of Helmi et al. (2006), namely that the MW 
halo was not formed from sub-haloes similar to the cur- 
rently observed satellites, may not hold either. Indeed, if 
an early infall phase is assumed for all the constituent 
sub-haloes, the composite halo DMD will necesarily dis- 
play a deficit of low metallicity stars, in agreement with 
the observations. However, it is impossible at present to 
determine the characteristics of the infall for each one of 
the sub-haloes as a function of their mass. 

For illustration purposes, we adopt here a simple toy 
model, where each sub-halo undergoes an early infall 
phase (of primordial composition), with the smaller sub- 
haloes having a more important infall. This prescription 
for infall is motivated from the observed DMDs of the 
satellite dSphs of the MW: as can be seen in Fig. 2, the 
most massive satellites, like e.g. Fornax, are better de- 
scribed by the simple outflow model, whereas in the low 
mass satellites the simple model fails badly in the low 
metallicity tail and requires an important early infall. Of 
course, in all cases, outflow at an appropriate rate (i.e. as 
to reproduce the observed effective yield) is assumed. 

The results of the toy model appear in Fig. 6. Because 
of the assumed early infall, the low metallicity tails of the 
small sub-haloes are less pronounced than in the simple 
model. Since those sub-haloes dominate the local DMD 
(because of their large number and low metallicity peaks) , 
the low metallicity tail of the composite DMD is also af- 
fected: the corresponding curve flts now the (preliminary) 
Hamburg/ESO survey data. Thus, it is possible that all 
types of galaxian systems, namely both the components 
of the MW halo and the small dSphs of the local group, 
have undergone an early period of infall, shaping the low 
metallicity tails of their DMDs. 

4. Summary 

In this work we present an analytical model for the metal- 
licity distribution of the MW halo, in the framework of 
the hierarchical merging scenario for galaxy formation. It 
is assumed that the MW halo was composed of a num- 
ber of sub-haloes with stellar masses ranging from Mi =2 
10^ Mq to M2=2 10* Mq and distributed as in Eq. (9). 
Furthermore, each of the stellar sub-haloes has a DMD 
given by the simple model (Eq. 2) with a mass depen- 
dent effective yield given by Eq. (4). Those assumptions 
are based on observations of dSph satellites of the MW 
galaxy, while the distribution function is obtained from 
numerical simulations of dark matter halo formation (af- 
ter explicitly accounting for modiflcation due to outflows). 

It is found that, under those assumptions, the bulk of 
the MW halo DMD is reproduced in a fairly satisfactory 
way. The advantage of that model over the simple outflow 



N. Prantzos: HMD 



9 



' I ' ' ' ' I ' ' ' ' I ' ' ' ' I ' ' ' ' I ' ' ' ' I ' ' 




-5 -4 -3 -2 -1 
[Fe/H] 



Fig. 6. Differential mctallicity distribution of the compos- 
ite MW halo population. The difference with Fig. 5 is that 
an early infall phase is assumed for each one of the com- 
ponent sub-haloes (except the most massive ones, see Sec. 
3.5), as to produce a paucity of low mctallicity stars. The 
dotted curve is the simple outflow model; data is as in Fig. 
1. 

model (proposed by Hartwick in 1976) is that it treats the 
MW halo not as a single entity but as a composite one, 
in line with the hierarchical merging scenario for galaxy 
formation. Furthermore, the proposed composite model 
avoids the "paradox" noticed in Sec. 3.2, namely that the 
effective yield of the MW halo appears to be substantially 
lower than the effective yields of dwarf galaxies of compa- 
rable mass. 

Observations suggest that the low mctallicity tail 
of both the MW halo and local dSphs is defficicnt in 
stars w.r.t. the predictions of the simple outflow model. 
Prantzos (2003) suggested that early infall may account 
for that defficiency in the case of the MW halo, while 
Helmi et al. (2006) prefer pre-enrichment (at a level of 
[Fe/H]=~-2.8) for the dSphs. It is shown here, with a 
simple toy model based on observations of dSphs, that 
an early infall phase can account for the low mctallicity 
tails of all small galaxian systems and, consequently, of 
the MW halo conceived as a composite system. 

The results obtained here with semi-analytical mod- 
els should be substantiated with detailed numerical simu- 
lations of the formation of a stellar halo in a MW-type 
galaxy. Such simulations suffer at present from several 
drawbacks, which make it difficult for them to provide 
robust predictions. For instance, the issue of star forma- 
tion vs feedback is far from being settled. Furthermore, as 
shown (and stressed) in Prantzos (2003), accounting for 
the finite lifetimes of stars produces substantially different 



results than the commonly adopted Instantaneous recy- 
cling approximation (IRA) at the earliest epochs (and low- 
est metallicities) ; since almost all numerical models adopt 
IRA (exceptions being Samland et al. 1997 or Lia et al. 
2002), their predictions cannot be considered as accurate, 
even if the issue of feedback is satisfactorily resolved. 

Finally, as emphasized in Sec. 2, an in-depth discussion 
of the important implications of the low mctallicity tail of 
the DMD should rely on the final release of the results of 
the Hamburg/ESO survey (Schorck et al. in preparation); 
the toy model presented here relies on preliminary results 
and should only illustrate the potential of the proposed 
analytical model to probe the physics of the early galaxy. 
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